Bound states in a nonlinear Kronig Penney model 



Stavros Theodorakis and Epameinondas Leontidis 

Dept. of Natural Sciences, University of Cyprus 
P.O. Box 537, Nicosia 1678, Cyprus 

February 1, 2008 



Abstract 

We study the bound states of a Kronig Penney potential for a nonlinear 
one-dimensional Schrodinger equation. This potential consists of a large, 
but not necessarily infinite, number of equidistant (^-function wells. We 
show that the ground state can be highly degenerate. Under certain con- 
ditions furthermore, even the bound state that would be normally the 
highest can have almost the same energy as the ground state. This holds 
for other simple periodic potentials as well. 



1 Introduction 



In this paper we shall study an unusual generalization of the one-dimensional 
Kronig Penney model. We shall examine in particular the spectrum of the bound 
states for a Kronig Penney potential V(z), having added though a nonlinear term 
to the Schrodinger equation. Our arguments will be valid in the case of other 
simple periodic potentials as well. 

Such nonlinear equations with periodic potentials arise in the Ginzburg-Landau 
treatment of various phenomena in condensed matter physics. In layered super- 
conductors for example, such as the high temperature ones, a periodic potential 
such as the Kronig Penney potential can describe the periodically modulated 
superconductivity of the samples [1]. Spatially varying parameters in the non- 
linear Schrodinger equation were also used to describe the periodic variation of 
the impurity concentration in superconductors [2], high T c Josephson field effect 
transistors [3], as well as grain boundaries in superconducting bicrystals [4], while 
nonlinear Kronig Penney models were used for studying twinning-plane super- 
conductivity [5]. The nonlinear Schrodinger equation must be used in order to 
describe all these various phenomena, including the relevant phase transitions. 
The nonlinear Schrodinger equation has been studied repeatedly, but mostly in 
regard to its solitons [6] , and usually for nonperiodic potentials. In this work the 
emphasis is placed on studying the bound states, rather than solitons. 

We shall study the excited states for the equation 

h 2 d 2 v 



+ V(z)ty + f3\^\ 2 ty = 0. (1) 



2M dz 2 

The nonlinear term forbids the arbitrary normalization of 

The potential we have in mind is a Kronig-Penney potential, but it could be 
in general any simple oscillatory potential. In this work we choose 

V(z) = V [l-aj:5( Z --n- 1 -)}, (2) 

with a and Vq positive. The crucial parameter in this potential is the periodicity 



length d. The number of wells is large, but not necessarily infinite. 
Equation (1.1) minimizes the energy functional 

J dz\v(z)\<S>\ 2 + p\^\ 4 /2 + 
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For M — » oo we would have |\I/| 2 = —V(z)//3, in which case |\l/| 2 would follow the 
periodicity of V(z). If the nonlinear term is omitted, the usual Kronig Penney 
model is recovered. In that limit — Vq is the energy, and aVo is the strength of 
each attractive delta function. 

We can write the energy functional in dimensionless form, by measuring z 
in units of d, the distance between successive spikes of the potential, VP in units 



of \JV /{3, and the energy in units of dV^/ft, where V is the positive constant 
that appears in Eq. (1.2), and has the dimensions of V(z). This constant is 
taken out of V(z), so as to render it dimensionless. In other words, V(z)/Vq = 
u(z), where u(z) is dimensionless. If we define then the dimensionless parameter 
v = h 2 /2MV d 2 , the energy functional takes the dimensionless form 

dz\u(z)\*\ 2 + m*/2 + v& 2 ]. (4) 
L dz J 
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Note that when the quartic term is omitted, we recover the usual linear Kronig 
Penney model, with energy E=-h 2 /2Mud 2 . In this case the energy values can 
be found only after imposing periodic boundary conditions on l^] 2 . There are 
then only certain allowed values of u, for a given value of a. The size of the 
wavefunction is determined by the normalization, and when we minimize the 
energy functional under this constraint, we find the energy eigenvalues, i.e. the 
minima of the energy functional. 

In the nonlinear case on the other hand, the size of ^ is determined by the 
nonlinear terms, through the unconstrained minimization of the functional of Eq. 
(1.4). These nonlinear terms determine fully the behaviour of VP, without any 
need for boundary conditions. In fact, a periodic u(z) will give a periodic |\P| 2 . 
Furthermore, the parameters a and v are now independent, and for any pair of 



values of a and v we can find a solution as long as a is sufficiently large. We 
shall see later what is precisely the lower bound on a. The energy of each state 
will be simply the value of the energy functional (1.4) at its minimum. 

We see from Eq. (1.4) that for v — > 0, when the potential is very strong, 
or very weakly periodic, we get l^l 2 — >— u(z). Thus \& follows very closely the 
periodicity of the structure, since it can change very abruptly. In this limit though 
the sign of \I/ is arbitrary. So if the spikes of the potential are very far apart (d 
is long), the sign of ^ could be positive or negative at each spike (see Fig. la). 

Let us now switch on slowly the parameter u, bringing the teeth of the po- 
tential comb closer together. Then the wavefunction between neighboring spikes 
could have two forms. If the wavefunction on two successive spikes A and B is 
positive, say, then the wavefunction in the intervening region will be reduced, and 
it will go through a positive minimum value, remaining always positive though 
(Fig. lb). If however the wavefunction changes sign in going from spike B to 
spike C, then it must pass through a point halfway between the spikes where it is 
exactly zero (see Fig. lb). Since the wavefunction ^(z) minimizes the functional 
of Eq. (1.4), the energy equals — / dz |\E'| 4 /2, as can be deduced by combining the 
dimensionless forms of Equations (1.1) and (1.3). Consequently the wavefunction 
has less energy if it does not go through zero, maintaining always the same sign. 
Indeed, in that case the minimum of |\l/| 4 /2 is not zero, and hence the area under 
|\E'| 4 /2 is greater. 

It seems therefore more favorable for the wavefunction to have the same sign 
on all spikes of the potential. We say that the ground state is a uniformly 
positive state then. If the spikes of the potential are too far from each other 
however, then the minimum value of the wavefunction between them is practically 
zero, and in that case the uniformly positive state (where \& has the same sign 
at all spikes) becomes degenerate in energy with states that may have \I/ take on 
negative values at some spikes, and positive values at others. 

We can, for example, have a state that is infinitesimally higher in energy 



compared to the uniformly positive ground state, and hence practically equally 
preferable, even for spikes not too far apart. This state, with \l/(0) = 0, ^(z) > 
when z > 0, and ^(z) < when z < 0, connects regions of different signs of the 
wavefunction (see Fig. 2). Then in the intermediate region \l/ has to go through 
zero, and we get a region that reminds us of a domain wall. For a potential with 
infinitely many spikes, the energetically costly root of ^ occurs only once, and 
hence the energy of this state is equal to the energy of the uniformly positive 
ground state. 

There can be bound states of Equation (1.1) therefore that are degenerate to 
the ground state, not being everywhere positive. It is the purpose of this paper 
to study such bound states, first through a general variational model (Section 
2), and then through an exact study of the Kronig- Penney potential (Section 3), 
as well as through a numerical study of a periodic potential with gaussian wells 
(Section 4). We summarize our conclusions in Section 5. 

2 Variational Study 

In this section we shall examine the possibility of having \I/ change sign in going 
from one spike of the potential to the next, as well as the possibility of having 
\& keep the same sign on neighboring spikes. In the first case \I/ is odd with 
respect to the midpoint between the two spikes, while in the second case it is 
even. An arbitrary state of the system will then be a combination of even and 
odd pieces. In other words, \I/ will be even between certain neighboring spikes of 
the potential, and odd between others. Thus ^ will maintain its sign between 
some spikes, and it will change sign between others. For example, in Fig. lb 
\l/ is odd in one interval, and even in the other two, while in Fig. 2 it is even 
everywhere, except for the interval at the center. 

We have assumed that the spikes of the potential are at z — n + |, where 
n is any integer. Let us examine then the two neighboring quantum wells at 
the ends of the interval [n — \, n + |]. We adopt the following odd and even 



trial wavefunctions, with respect to the midpoint (z — n), defined on the interval 
[n - ±,n+ \\: 

smhhiz - n)] 

(z) = ±*P n , 5 

smh(7/2) 



V en (z) = ±V> 



coshMz — n)] , o , 

J - sech 2 ( 7 /2) 



coth 2 ( 7 /2), (6) 



cosh(7/2) 

where tp and 7 are variational parameters. We note that \l/ ori (n) = 0, ^ on (n + |) 
= -^on(n-|) = ±tp, and % n {n + \) = % n (n-±) = ±^7 coth( 7 /2). Similarly 
V e n(n) = ±^[cosh( 7 /2) - l]/sinh 2 (7/2), ^ en (n + \) =^ en (n - \) = ±ijj, and 
Knin + \) =~Ka in-\) = ±^ 7 coth( 7 /2). 

These wavefunctions are such that they can be joined together in any order to 
form a continuous wavefunction everywhere, consisting of even and odd pieces. 
We could have, for example, ^ = \l/ on in [n— ^ = ^n+i in [n+|,n+|], 

^ = —^ ,n+2 in [n + |, n + |], etc. Furthermore, regardless of the order in which 
the even and odd pieces are connected, the slope of the wavefunction is symmetric 
around the spikes of the potential. 

The state with the lowest energy would consist of a chain of even pieces, 
because, unlike the odd pieces which have a root at the midpoint, the even pieces 
are nowhere equal to zero. Thus the odd pieces have a higher — / dz\^\ A /2, 
which is the exact energy if \1> is an exact solution of the equations that minimize 
the energy functional. 

The next lowest energy would correspond to the state with only one odd piece. 
This is the state of Fig. 2. It is presumed here though that the change from the 
chain of negative pieces to the chain of positive pieces occurs within just one 
spacing. The circumstances under which this will happen will be explored later. 

The state mentioned above is followed by the state with two odd pieces, and 
so on, up to the highest state, which has only odd pieces. In fact, if F e and F G is 
the energy in [n — |, n + |] for the even and odd trial wavefunctions respectively, 
then the energy of a state with m even pieces and n odd pieces is mF e +nF D . Thus 
the total energy per interval is (mF e +nF )/(m+n). In particular, the energy per 



interval is F e for the uniform chain of even pieces, i.e. the ground state, and 
(mF e +F )/(m+l) for the state with only one odd piece, i.e. the state of Fig. 
2. For an infinite number of spikes (m— > oo), the two states are degenerate, 
as expected. Of course, the same holds for a state with infinitely many even 
pieces, but only two odd pieces. If the number of odd pieces becomes substantial 
though, then the energy of the state will be definitely higher than that of the 
ground state. 

For a large but finite number of spikes we still expect all these various states to 
be degenerate, as long as the minimum value of the even pieces is practically zero, 
because in that case |\1/| 4 is essentially the same for both even and odd pieces. 
We shall verify this by explicit calculation, using our variational wavef unctions. 

We note that ty en (n) — > when 7 — > 00. Therefore the degeneracy men- 
tioned above requires that 7 be very large. We shall neglect therefore terms like 
sech(7/2). In this limit, 

r.~rW + & + £! ^MZg^ + ow (S) 

Hence, since cosh 2 (7/2) sinh 2 (7/2) ^e 7 /4, 

F -F e ^- dz4e- y u(z)\ij\ 2 . (9) 

Jn-\ 

And since 7 is large, and terms of order 0(e -7 ) have been dropped in this cal- 
culation, Eq. (2.5) implies that F D ~F e . In other words, if \& en (n) ~ 0, then all 
the possible states are practically degenerate, even for a finite large number of 
spikes, because they consist of odd and even pieces only, pieces which were shown 
to have the same energy. Note that our results are very general so far. The only 
restriction is that the 7 that minimizes Eq. (2.3) be large. Our conclusions are 
valid though for any u(z) that can lead to a large 7. 

We illustrate the above general conclusions by restricting ourselves now to the 



Kronig-Penney model: 



u{z) = l-^2aS{z-n--). 



(10) 



This choice of u(z), where a is a positive constant, implies that there is a periodic 
chain of deep quantum wells along the z axis. 



For this choice of u(z) then, and in the limit of large 7, we get 
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Minimization with respect to j^] 2 gives 

M 2 = 2 7 
and 



a — 1/7 

7- 



F e « -7 



ct — 1/7 — 



7 



Minimization with respect to 7 yields 

a + \Jo? + 12z/ 



7 



61/ 



(11) 
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(13) 



(14) 



or equivalently 3z/7 — a = I/7, in which case \ip\ 2 = Awy 2 — 4. Since 7 ^> 1, we 
shall have 

7«q/3i/»1. (15) 
Since j-?/'! 2 > 0, we must also have 7 > l/y/v, which implies 



a 2 > Au. 



(16) 



Whenever therefore v and a satisfy the restrictions of Eqs. (2.11) and (2.12), 
we expect all the possible states to be essentially degenerate. In particular, the 
highest excited state, the one consisting of odd pieces only, is degenerate with 
the ground state, which is a chain of even pieces. Note furthermore that the 
wavefunction is nonzero only if v < a 2 / 4. When in fact v = a 2 / '4, we have 
a transition to a zero wavefunction, even though 7, which takes then the value 
1 / y/u, may be quite large. 

As a numerical illustration, we choose the case v = 0.01, a — 1. Then 
the ground state and the highest excited state (only even or only odd pieces 



respectively) have an energy of —13.5207 in this variational model, with 7 = 
34.305 and ip = 6.56. The exact energy can be found using the methods of 
section 3, and it is —14.933 for both the ground state and the highest state, 
while ip = 6.93. So both calculations indicate that all the states are degenerate, 
for this particular choice of v and a. 

The example where u(z) is given by Eq. (2.6) will be examined also in section 
3, since it can be solved exactly. We can generalize for any simple oscillatory 
u(z) though, so long as 7 is very large. 

3 Exact Solutions 

In this section we shall solve exactly the model of Eqs. (1.4) and (2.6), verifying 
thus the variational results of the previous section. We shall be interested in those 
values of the parameters v and a that yield excited states almost degenerate with 
the ground state. 

We should note that a large v would imply that the kinetic energy is domi- 
nant, making thus the wavefunction too stiff. In other words, the ground state 
wavefunction would come as closely as possible to a constant, a choice that min- 
imizes the kinetic energy. In that case the minimum value of the ground state 
wavefunction would be far from zero. 

On the other hand, if v were exactly zero, then the wavefunction would follow 
the variations of u(z) exactly. Hence we need a small value of v if we are going 
to have an excited state that is close in energy to the ground state, since the 
wavefunction of such a state varies dramatically between the spikes. Furthermore, 
if v is zero the wavefunction will have arbitrary signs at the wells, in which case 
the various excited states will all be degenerate with the ground state. For small 
is, this degeneracy will not be altered too drastically. 

We shall be interested therefore in the exact solutions of this model, for small 



v. The energy functional is minimized when 

d 2 m 



v 



dz 2 



(17) 



The solution \I/(-2) will have periodic features similar to those of u(z). Integrating 
Eq. (3.1) gives the boundary condition for ^f(z) 

" + 2) = A T z ^ + 2 )+ - + 2 } -l • (18) 

Thus *&(z) has a kink at each spike of the potential, due to the ^-functions. 

Direct integration of Eq. (3.1) after multiplying it by /dz gives the solution 
in each interval. The ground state has no node, hence \& will have a minimum at 
the middle of each interval, while it will be symmetric around each spike. Thus 
the exact ground state is found to be 



Q 



cn 



(19) 



[l + q 2 )/u (z-n), (2 + g 2 )/(2 + 2g 2 ) 

for n — I < z < n + |, extended periodically everywhere else. Here cn is a 
Jacobi elliptic function, and q = ty(n) is the minimum value of ^(z). The above 
expression is valid for any value of is, large or small, and we can easily verify that 
it satisfies Eq. (3.1). 



The boundary conditions of Eq. (3.2) require then that 



a = 2u\ 



h + q 2 sn 


\/lW 2W~ 
L V Av ' 2+2q 2 _ 


dn 


_V Av ' 2+2y 2 _ 




cn 


/1+q 2 2+q 2 ' 

_V Av i 2+2q 2 _ 





(20) 



where the dn and sn are also Jacobi elliptic functions. This equation determines 
q as a function of a and v. Note that there is always a ground state, since we 
can always find an appropriate q for a given choice of v and a. 

In the limit of a wavefunction localized around the spikes of the potential we 
expect q to be small. Then Eq. (3.3) reduces to 



m gcosh[(,2 — n)/y/u\. 



(21) 



This is precisely the solution of the linear Kronig Penney model, as expected, 
since for very small q the nonlinear terms become unimportant. If q — 0, Eq. 
(3.4) yields 

a = 2v^tanh(l/2v^). (22) 

So a nonzero ground state will exist only if a > 2 v /z/tanh(l/2 v /z/). Note that 
for small v this becomes the restriction of Eq. (2.12), as it should. Furthermore, 
if a is close to its lower limit, then the change of sign of ^ for a first excited 
state will have to spread over a few more spikes of the potential, and it will not 
be restricted to just the region between two successive spikes. In this paper we 
shall not be concerned with this possibility, and we shall restrict our attention to 
values of a far from the lower bound of Eq. (3.6). Then the change of sign for 
the excited states occurs within just one spacing. 

Now the Jacobi elliptic function cn(x, m) is a periodic function, with roots at 
the odd multiples of the elliptic function K(m), where K(m) = f^ 2 d&j\J\ — m sin 2 9. 
Indeed, cn(0, m) = 1, cn(K(m),m) = 0, cn(2K(m),m) = — 1, cn(3K(m),m) = 
0, cn(AK(m),m) = 1. 

The ground state, as mentioned in the previous section, must consist of even 
pieces everywhere, of the form given by Eq. (3.3). It must therefore be positive 
everywhere, since it is continuous. So the \& of Eq. (3.3) must not be allowed 
to become negative. This means that the quantity \J(1 + q 2 )/v\z — n\ must be 
smaller than K{{2 + q 2 )/(2 + 2q 2 )) within the interval [n — \,n + |]. Hence 

\y/t&W»<K{££). (23) 

This inequality holds for any values of a and v. If v is small, then the right hand 
side of this inequality has to be large. This happens when the argument of K{m) 
is close to 1, in which case K{m) In ^16/(1 — m). In this particular case, the 
argument is 1 if q is very small. Thus if v is small, q must be small. 

Furthermore, when inequality (3.7) becomes an equality, the denominator in 
Eq. (3.3) tends to zero at z — n + 1, and therefore the value of *&(z) at the spikes 



becomes infinite. In other words, the value of \1/ at the spikes can be much larger 
than the value of \I/ at the midpoints. 

Indeed V(n - (n) = l/cn[^J(l + q 2 )/4u, (2 + q 2 )/(2 + 2q 2 )]. Hence, if 
q — > 0, ^(l/2)/^(l) ks cosh(l/v^4z/), which tends to infinity when v tends to 
zero. Thus the even pieces of ^ become very deep if v is small, because in that 
case q — > and \l/(n ± |) is very large. 

We say in that case that the wells are weakly coupled. Note that in that case 
the |\l/| 4 of an even piece would not differ too much from the |\l/| 4 of an odd piece. 
In other words, we expect the various possible states to be very close in energy 
to the ground state, as mentioned already in section 2. Indeed the case v — 
would correspond to a complete decoupling of the values of the wavefunctions at 
the spikes of the potential, and hence to a complete degeneracy of all the various 
states. 

Let us examine more thoroughly the singularities that may arise in the be- 
havior of We said that becomes very large at the spikes of the potential 
when inequality (3.7) becomes almost an equality: 

y/(l + ^)/Au^K(^±^). (24) 

But if m is very close to 1, then K(m) « In ^/l6/ (1 — m). Hence this approximate 
equality reduces for small q to l/2y/v ~ ln(y / 32/g 2 ), whence 

q « Vzle- 1 ' 2 ^. (25) 

Hence, if v is small, and if the ground state wavefunction has deep cups, we must 
have q xs \/32 exp(— l/2y/v). Note that even though q is small, the value of the 
wavefunction at the minima of the potential is large. Nonetheless, since q — > 0, 
the wavefunction is again given by Eq. (3.5), an equation that tells us that ^(z) 
falls to 1/e of its value within a distance of \Jv from the spikes. In that sense we 
can say that the "thickness" of ^ at each spike is But the peaks of \P would 
overlap when the thickness of each peak equals the distance between successive 
peaks. This happens when v — 1/4. When we speak therefore of weakly coupled 



wells, we mean that v <C 1/4. And it is only such wells that can lead to an 
essentially degenerate spectrum of states. 

Let us then summarize our results for the ground state. There is always 
a ground state, with q = ^(n) being the minimum value of \1/ in the interval 
[n — |, n + |]. This is a symmetric series of even pieces, and it resembles a chain 
of symmetric cups (see Fig. 3). The absolute value of the slope of \1/ at the layers 
is a^/(n + |)/2za This ground state will be degenerate with any other states 
if their |\l/| 4 's are approximately the same. This can happen only if q is almost 
zero, as explained in section 2, because then the minimum of the |\1/| 4 of the even 
piece approaches the minimum of |\&| 4 of the odd piece, i.e. zero. But q can be 
tiny, and \l/ still have a substantially nonzero value, only close to the roots of the 
Jacobi elliptic finction cn (see Eq. (3.3)), i.e. for q « v^e" 1 / 2 ^ 7 . Furthermore, 
q needs to be small in order to have the degeneracy. Hence v must be small. 
Indeed, the thickness of each well implies that the wavefunctions around the 
spikes will not overlap substantially, provided v <C 1/4. For small v then we get 
a ground state which resembles a chain of deep cups (see Fig. 3) . 

Let us now proceed to the first excited state (see Fig. 2). Here we assume 
again that v is small, and hence the even pieces will resemble deep cups. There 
will be only one odd piece, in the interval [—1/2,1/2], connecting a chain of 
negative even pieces with a chain of positive even pieces. The characteristics of 
the many even pieces will not be altered, because there is only one odd piece. On 
the contrary, the characteristics of the odd piece will be determined from those 
of the even pieces, through the boundary conditions. 

Direct integration of Eq. (3.1) after multiplying it by d^/dz will give the 
solution in the interval [—1/2, 1/2], as long as we use the fact that ^(O) = 0, 
since there is one node there. For the first excited state there is only one node, 



thus all the pieces outside the interval [—1/2, 1/2] will be even. 

One can show thus that for —1/2 < z < 1/2 we get the exact solution 



, sn[zJ(l + e)/2u, 2e/(l + e)] 

y( z ) = JY^l v ; =— -, (26) 

cn[zJ(l + e)/2i/, 2e/(l + e)] 



with #'(0) = ^/(l - e 2 )/2z/ and < e < 1. We can easily verify that this 
expression satisfies Eq. (3.1). The value of #(2) at z — 1/2, as calculated from 
Eq. (3.10), must be equal to the one that can be calculated from Eq. (3.3). This 
relation determines the parameter e. If the even pieces are deep enough, i.e. if v 
is small enough, then the slope of ^(z) at (l/2)_ will turn out to be a^{\)/2v. 

Indeed, we saw that if for small v the value of # at the spikes of the potential is 
very large, then q ps v^e -1 / 2 ^. In general, the even and odd pieces correspond 
to the same energy if # is very large at the spikes, q being quite small. But if 



#(1/2) is very large, then Eq. (3.10) implies that cn[yj{\ + e)/8is, 2e/(l + e)] « 0, 
so as to make *&(z) almost diverge. Consequently 

y/(l + e)/8vttK[2e/(l + e)], (27) 

where K is the elliptic function K(m) = J^ 2 d9/\/l — msin 2 6. And since v is 
small, the value of K will have to be rather large, which means that 2e/ (1+e) — > 1, 



i.e. e — > 1. Indeed, remembering that K(m) ~ In y 16/ (1 — m) when m — > 1, we 
can easily find that Eq. (3.11) is solved by the value 

1 - e w 32e~ 1/v ^. (28) 

So this value of e yields a very large #(1/2), for small v. In fact, we must have 
in general, for any u, 



;i + e)/8i/< tf[2e/(l + e)], (29) 

otherwise the elliptic function cn would get a root in [0,1/2] and #(z) would have 
a vertical asymptote there. 

We can now check the value of #(1/2). The second argument of the elliptic 
functions sn, dn and cn is (2 + q 2 ) /(2 + 2q 2 ) for the even pieces and 2e/ (1+e) for 



the odd pieces, both of which will equal 1 — lQe^ 1 ^ when Eqs. (3.9) and (3.12) 
hold, i.e. for very large values of \I/(l/2). Therefore for very small v this second 
argument is essentially 1, in which case the sn becomes tanh, the cn becomes 
sech, and the dn sech. Then Eq. (3.10) gives *&(z) « y/1 — e smh(z\J (1 + e)/2v) 
in [—1/2, 1/2], and qcosh\(z — n)/y/v\ in [n— \,n + \], where n^O. Both 

expressions give then the same values for ^(1/2) and |^'(1/2)|, as expected. 

The procedure for finding the first excited state then consists of finding the 
value of e that would ensure continuity of ^(z) at z = 1/2. In that case the slope 
at z — (l/2)_ will turn out automatically to be the exact opposite of the slope at 
z = (1/2)+. Finding the ground state, on the other hand, simply requires finding 
a q such that |\&'(n + |)| = (a/2v)^( y n + |). This relation is precisely Eq. (3.4). 

Finally, we can find the highest excited state, the one consisting of odd pieces 
only, by extending periodically the odd solution of Eq. (3.10), and finding a value 
of e such that \^'(n + |)| = (a/2v)^(n + |). Since the minimum value q of the 
even piece for the solutions that interest us is v^e^ 1 / 2 ^, i.e. practically zero, 
the energy of the even piece and of the odd piece is essentially the same since 
they have the same — / dz\^\ 4 /2. And all the states are then degenerate. 

It is interesting to note that a = 2i/\l/ / (i)/\l/(|)— > Av + | if v — > oo, for the 
^(z) of Eq. (3.10). Hence the highest state does not exist when v — > oo, unless 
a = 4z/ + |. Similarly, quite a few other excited states do not exist for large 
values of v. The even ground state exists always though, if a is above the lower 
bound of Eq. (3.6). We present numerical values of the parameters and energies 
of the ground state, of the first excited state and of the highest state in Tables 
1, 2 and 3 for various choices of v and a. 

Note that for a given choice of v, large or small, and a given value of a greater 
than the lower bound of Eq. (3.6), q is given by Eq. (3.4). But it is only small 
z/'s, and a's quite far from their lower bound, that will lead to a ground state 
with deep even pieces. In that case q is given by Eq. (3.9), and the ground state 
will be practically degenerate with the first few excited states. 



4 A Potential With Gaussian Wells 



We shall examine numerically the bound states that correspond to the choice 

u(z)= 1 - g^eM-Kz-n- 1 -) 2 }, (30) 

n Z 

shown in Fig. 4. We choose the parameters so that e fe//4 > g > 1. Then u(z) 
is negative at its minima and positive at the midpoints between the minima. 
Therefore, if v were exactly zero, then there would be a nonzero wavefunction at 
the minima of u(z), but there would be a region around the midpoints where \I/ 
would be exactly zero. When v is slightly positive, the stiffness of the wavefunc- 
tion makes \l/ leak into the "forbidden" regions (tunnelling). If however b < 41ng, 
then there is no tunnelling. We need a large enough b in order to get tunnelling. 

If on the other hand b is extremely large, while g remains finite, the width of 
each well is reduced, its strength remaining unaltered. We expect therefore that 
for a given g we cannot increase b indefinitely, because we shall not be able to 
find a solution. The numerical calculations do indeed verify this. 

Thus, if b is large enough, but not too large, and if g is large enough to allow 
a nonzero solution, we shall have a competition between the kinetic energy and 
the rest of the energy. The latter forces \1/ to follow the variations of u(z). The 
kinetic energy, on the other hand, wants ^ to be constant in space. Solutions 
such as the one of Fig. 2 will be relevant for small z/, since large values of v 
tend to push ^ towards a constant. Indeed, numerical calculations verify that 
for large v (e.g. v = 0.3), the kinetic energy is strong enough to force ^ ma x to be 
close to ^min, and ^ m i n ^> 0. The ground state is then a chain of shallow little 
cups away from zero. 

Excited states that will be almost degenerate in energy with the ground state 
will appear at small values of u, when the minimum of the ground state wave- 
function between the wells is very small, while $ mm ^> 0. Fig. 5 shows the 
ground state, the first excited state, and the highest excited state, for v = 0.01, 
with the corresponding energies per interval, for a sample of 20 gaussian wells. 



Fig. 6 shows the same states, but for v = 0.05. We see again the characteristics 
mentioned earlier: full degeneracy if v is small enough, in which case the mini- 
mum of \& between wells is quite small compared to its value ^ max at the wells. 
The value of ^ max depends strongly on g. 

Remember also that the energy difference between the ground state and the 
first excited state is strictly equal to zero for an infinite number of wells. This is 
due to the fact that the wavefunction of the first excited state is exactly equal 
or exactly opposite to the wavefunction of the ground state on almost all of the 
infinitely many wells. Indeed, in Figs. 7 and 8 we see the first excited state for a 
fairly large value of u, and a series of values for b and g. Here the stiffness of the 
wavefunction is quite large, so the change of the sign cannot take place within 
just one interval. This change now occurs over three or four wells. For small 
values of v though, when the wavefunction is quite malleable, the wavefunction 
changes sign within one interval. 

5 Conclusions 

We have studied the bound states of a nonlinear version of the Schrodinger equa- 
tion for the Kronig Penney model, a version that is relevant to quite a few phe- 
nomena in condensed matter physics. 

We have seen that there is a substantial range of parameters, not just for the 
Kronig Penney model but for other simple oscillatory choices of u(z) as well, for 
which the various states are essentially degenerate in energy. This degeneracy 
requires that v be small enough, so as to allow the wavefunction to have a small, 
but non-negligible, value at the midpoints between the wells. At the same time 
the value of the wavefunction at the wells can be quite substantial. Then there 
is very little cost in having \I/ change sign in going from one well to the next. 
In fact, the energy differences are really small when f mm ^> ^ m in- This means 
that the lower excited states, which connect regions of positive \l/ with regions of 
negative become as favorable as the ground state. This degeneracy is exact in 



the limit of infinitely many quantum wells for these lower excited states. Even 
the highest excited state though, which consists of odd pieces only, has an energy 
very close to the energy of the ground state for sufficiently small v. 

We demonstrated this basic idea in Section 2 through a variational calculation 
valid for a generic oscillatory potential. The variational results were confirmed 
through the exact solution of the nonlinear Kronig Penney model, presented in 
Section 3, as well as through a numerical calculation for the case of another 
simple oscillatory potential, presented in Section 4. All these calculations show 
the exact (for infinitely many spikes) or approximate (for finitely many spikes) 
degeneracy of the ground state. 
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Table Captions 



Table 1: Parameters of the ground state for 5-function wells, for various choices of v 
and a. E g is the energy per interval. 



V 


a 


*(1) = 5 


tf(l/2) 




E 9 


0.002 


1 


7.2xl0~ 5 


15.748 


3937 


-81.4526 


0.01 


0.25 


0.0127 


1.0608 


13.2601 


-0.02710 


0.01 


1 


0.031074 


6.92821 


346.41 


-14.9334 


0.01 


10 


0.037283 


70.6965 


35348.3 


-16646.9 


0.05 


1 


0.34647 


2.83402 


28.3402 


-1.95903 


0.05 


10 


0.511176 


31.5911 


3159.11 


-3314 


0.1 


1 


0.497576 


1.78192 


8.90962 


-0.59964 



Table 2: Parameters of the first excited state for ^-function wells, for various choices 
of v and a. E 1 is the energy in [—1/2,1/2]. For the other intervals, the 
energy is the E g given in Table 1. 



V 


a 


Q 


e 


¥(1/2) 


¥'(1/2)- 


E l 


0.002 


1 


7.2xl0~ 5 


l-5.2xl0~ 9 


15.748 


3937 


-81.4526 


0.01 


0.25 


0.0127 


0.999838 


1.0608 


13.2613 


-0.027094 


0.01 


1 


0.031074 


0.99903 


6.9282 


346.41 


-14.9333 


0.01 


10 


0.037283 


0.998598 


70.696 


35348.3 


-16646.9 


0.05 


1 


0.34647 


0.82356 


2.8340 


28.4417 


-1.92445 


0.05 


10 


0.511176 


0.363395 


31.591 


3159.12 


-3313.85 


0.1 


1 


0.497576 


0.424940 


1.7819 


9.28766 


-0.526821 



Table 3: Parameters of the highest state for 5-function wells, for various choices of 
v and a. Eh is the energy per interval. 



V 


a 


e 


#(1/2) 


#'(l/2)_ 


E h 


0.002 


1 


l-5.2xl0~ 9 


15.748 


3937 


-81.4526 


0.01 


0.25 


0.999838 


1.0605 


13.2565 


-0.02707 


0.01 


1 


0.99903 


6.9282 


346.41 


-14.9332 


0.01 


10 


0.998598 


70.696 


35348.3 


-16646.9 


0.05 


1 


0.82442 


2.8213 


28.213 


-1.89582 


0.05 


10 


0.363395 


31.591 


3159.11 


-3313.84 


0.1 


1 


0.516054 


1.6527 


8.26345 


-0.40157 


1 


4.34 


0.730536 


0.2518 


0.54642 


-0.00038 



Figure Captions 

Figure 1: (a) A typical wavefunction when the potential wells are very far apart [y xs 
0). 

(b) A typical wavefunction when the potential wells are closer together 

Figure 2: First excited state with the root at z = 0, for 5-function wells. 

Figure 3: The ground state, for 5-function wells. 

Figure 4: The potential u(z) for gaussian wells (g = 2.5, b = 20). 

Figure 5: The ground state, the first excited state, and the highest state, for 20 
gaussian wells with v = 0.01, g=2.5, b=20, with the corresponding energies 
per interval. 

Figure 6: The ground state, the first excited state, and the highest state, for 20 
gaussian wells with v = 0.05, g=2.5, b=20, with the corresponding energies 
per interval. 

Figure 7: First excited state for gaussian wells with v — 0.1, b — 20, and g=2.2, 2.5 
and 5. 

Figure 8: First excited state for gaussian wells with v — 0.1, g — 2.5 and b =5, 10, 
20, and 25. 



